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ABSTRACT 

We have calculated self-consistent equilibria of a collisionless relativistic electron-positron gas in the 
vicinity of a magnetic X— point. For the considered conditions, pertinent to extra-galactic jets, we find 
that leptons are accelerated up to Lorentz factors Tq = kcBqLE 2 /mc 2 ^> 1, where Bo is the typical 
magnetic field strength, £ = Eq/Bo, with Eq the reconnection electric field, L is the length scale of 
the magnetic field, and k ss 12. The acceleration is due to the dominance of the electric field over the 
magnetic field in a region around the A— point. The distribution function of the accelerated leptons is 
found to be approximately dn/dj cx 7 _1 for 7 < Tq. The apparent distribution function may be steeper 

than 7 _1 due to the distribution of Tq values and/or the radiative losses. Self-consistent equilibria are 
found only for plasma inflow rates to the A— point less than a critical value. 

Subject headings: acceleration of particles: magnetic fields: MHD : plasmas: jets — galaxies 



1. INTRODUCTION 

The observed radiation of most large scale extragalactic 
jets is due to incoherent synchroton radiation. In a num- 
ber of sources, M87 and 3C 273 for instance, the radiation 
lifetime of the electrons (and possibly positrons) is much 
less than the transit time from the central source (Felten 
1968). Thus, there must be mechanisms for the electron 
"reacceleration." 

Several mechanisms have been proposed to account for 
the reacceleration of the electrons. These include Fermi 
mechanism of particle acceleration (Christiansen, Pachol- 
czyk, & Scott 1976), Fermi acceleration in shock waves 
(Krimsky 1977; Axford, Leer, & Skadron 1978; Bell 1977, 
1978; Blandford & Ostriker 1978), stochastic electric field 
acceleration (Eilek & Hughes 1990) and whistler wave ac- 
celeration (Melrose 1974). The effectiveness of the Fermi 
mechanism and stochastic fields in accelerating electrons 
is unknown (Blandford & Eichler 1987; Eilek & Hughes 
1990; Jones & Ellison 1991). Resistive tearing has been 
discussed as a means of producing neutral layers which can 
accelerate electrons (Koniglfe Choudhuri 1985; Choudhuri 
& Kdnigl 1986). Reconnection as a mechanism for accel- 
erating electrons has been discussed by a number of au- 
thors (Blandford 1983; Begelman, Blandford, k, Rees 1984; 
Browne 1985; Ferrari 1984, 1985; Norman 1985; Kirchner 
1988; Lesch 1991; Romanova & Lovelace 1992). The ac- 
celeration of electrons by long wavelength electromagnetic 
waves trapped in the boundary layer of a jet was discussed 
by Bisnovatyi-Kogan and Lovelace (1995). 

The reconnection of magnetic fields at a neutral point 
as described by magnetohydrodynamics (MHD) has been 
studied extensively. It has been shown that the neutral 
lines can evolve into current sheets (Syrovatskii 1971). 
Resistive MHD simulations of reconnection by Biskamp 
(1986) show a strong tendency to form a current sheets. 

Collisionless reconnection has been studied by following 
particle orbits, assuming nonrelativistic particles and pre- 
scribed fields (Galeev 1984; Zelenyi et al. 1984; Burkhart, 



Drake, & Chen 1990; Deeg 1991). The particles were 
launched away from the A— point and given an initial ve- 
locity that is appropriate for a given thermal velocity dis- 
tribution. Acceleration of the particles near the A— point 
was observed. This work was extended to include the self- 
consistent fields of the particle flows (Burkhart, Drake, & 
Chen 1991). For sufficiently high inflow rates they found 
an outflow shock where the flow velocity exceeded the 
Alfven velocity. In addition they found that sufficiently 
large inflow rates caused the width of the current layer to 
"collapsed to zero," and concluded that there was a maxi- 
mum infow rate into the reconnection region. Nonrelativis- 
tic particle trajectories have been studied in reconnection 
regions as part of an investigation of the Earth's magne- 
tosphere (Speiser 1965). This work has been extended in 
an effort to define an effective conductivity in the absence 
of particle collisions (Speiser 1970). 

Direct acceleration of particles in reconnection layers 
was proposed by Alfven (1968). Direct relativistic acceler- 
ation of leptons near the A— line of a magnetic field con- 
figuration was studied by Romanova and Lovelace (1992). 
The static electric field parallel to this A— line can cause 
the acceleration of leptons to the very high Lorentz fac- 
tors observed. Romanova and Lovelace (1992) analyzed 
the particle orbits in the current layer structure pro- 
posed by Syrovatskii (1971). For the case of an electron- 
positron plasma they obtained an energy distribution func- 
tion dn/dj cx 7~ 1,5 . Recently, Zenitani and Hoshino 
(2001) have done relativistic particlc-in-ccll simulations of 
a reconnection layer in an electron-positron plasma and 
have found a spectrum 7 _1 out to 7 ~ 20. 

In this work we solve the relativistic equations of mo- 
tions for an electron-positron plasma moving in the elec- 
tric and magnetic fields of a reconnection layer indicated 
in Figure 1. In turn, the particle motion is used to calcu- 
late the current-density and the associated self magnetic 
field. The particle orbits and field calculations are done 
iteratively to give a self-consistent equilibrium for the re- 
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Fig. 1. — Geometry of magnetic field reconnection layer. 



connection layer. Once we have this equilibrium we can 
determine the distribution function of the accelerated lep- 
tons. Section 2 of this paper develops the theory, and 
Section 3 discusses the methods and results. Section 4 
gives the conclusions of this work. 

2. THEORY 

2.1. Physical Picture 

The geometry of the stationary (d/dt = 0) reconnec- 
tion region is shown in Figure 1. The magnetic field 
B = B x (x,y)x + B y (x 1 y)y has an A— point in the (x, y) 
plane at x = 0, y = 0. A uniform electric field exists in the 
z-direction, E = E z with V x E = and E < 0. There- 
fore, the lepton gas drifts with velocity Vd y — cE z B x /Ti 2 
in the y— direction towards the magnetic A — point from 
above and below. Because the gas is electrically neutral 
there is no net current due to this drift (both electrons 
and positrons drift with the same velocity in the same di- 
rection). 

The B field vanishes at the A— point, and consequently 
in the vicinity of this point the particle motion is domi- 
nated by the electric field. The electric field accelerates 
electrons in say the +z direction and positrons in the — z 
direction. This gives a current-density in the — z direction. 
As the leptons are accelerated in the ±z directions their 
motion is "bent" into the ±x directions by the (v x B/c)^ 
force. At a large enough distance \x\ from the A— point the 
magnetic field becomes dominant and the particles exhibit 
the usual drifts (Northrop 1963). 

The electrons and positrons are confined near y — by 
the B x magnetic field. If a particle drifts into the X point 
from above, and then overshoots the y = plane, the 
E x B drift will then deaccelerate the particle and return 
it toward the y = plane. The particles also drift in the 
re— direction with velocity Vdx = —cE z B y /H 2 . Both elec- 
trons and positrons drift away from the A— point. This is 
due to the component of the B field in the y— direction. 
Since the a;— drifts for electrons and positrons are in the 
same direction there is no net current in the x— direction. 
Since both positrons and electrons enter the A— point from 
both sides (y > and y < 0) the net contribution of the 
y— motions cancels, and there is no y— direction current. 

We make the simplifying approximation of neglecting 
the y— thickness of the current layer. In this case, the 
particles move in the y = plane. We therefore treat 
equilibria where the particles drift from above and below 
into the y = plane where they are accelerated in the ±z 
directions and then expelled in the ±x directions. Thus 
we calculate the particle orbits [x(t), 0, z(t)\. Knowing the 
particle motion allows us to calculate the surface current- 
density J z {x). From this we calculate the self- magnetic 
field. We then use the self field to recalculate the par- 
ticle orbits. We iterate on this process so that we get a 
self-consistent solution of the self magnetic field and the 
particle orbits. The self-consistent orbits of the leptons 
can then be used to derive the energy spectrum of the 
accelerated leptons. 

2.2. Single Particle Motion 



The magnetic field can be written as B = V x 
p4 z (x, y)z], where A z is the total vector potential. For 
specificity we consider an A— type null point of the B so 
that A z (x,y) is an even function of both arguments. The 
magnetic field consists of a "external" component due to 
distant currents and the "self field" due to local currents. 
The external component A e z xt is divergence and curl free. 
We take the leading terms of a Tayor expansion of this 
field. 
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so that 

y_ 

'L ' ~ y ~"L 
where B > without loss of generality. The total field is 
given by 

A tot (x,y) = A ext + A self . (2) 

A tot is an even function of both arguments but its depen- 
dence on (x, y) is changed by the self-field. 

We consider quasi-stationary conditions so that V x E = 
and thus E = — V$, which is the "external" electric field. 
The relevant solution is E = and <& = —Eqz, where 
E = const < 0. This corresponds to plasma flowing into 
the A point from above and below the y = plane. Be- 
cause we consider an electron-positron plasma the "self 
electric field" is zero, and the total electric field is Eqz. 

The single particle motion is described by the La- 
grangian 



C = -mcV 1 - fi 1 + ~A z v z - q<S> 



where (3 = |v|/c and the electrostatic potential $ 
The corresponding Hamiltonian is 



(3) 
-E z. 



H 



(P - U) 2 + m 2 c 4 
c 



where P is the canonical momentum. Because dC/dt = 0, 
the H = const which is the single particle energy. 

A further constant of the motion is the canonical mo- 
mentum in the z— direction, P . Because A z is indepen- 
dent of z, 

f = "f -«*■ » >-^<> = o, 

so that 

dz Q 

P = P Z - qE Q t = m-y— + -A z - qE t = const . (4) 
at c 

The z— equation of motion is simply dP a /dt = which 
gives 

d ( dz\ q dx q dy „ „ ,_. 

sN + c B ^-^l-^ =o ' (5) 

where 7 = (1 — /3 2 ) -1 / 2 . For the limit where the 
z— acceleration is small, equation (5) gives 
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or 



v = cE x B/B 2 



(6) 



which is the well-known E x B drift velocity (Northrop 
1963). 

We assume that away from the X— point the leptons 
are relatively 'cold' in the sense that their thermal veloc- 
ity spread is less than their E x B drift velocity which is 
less than c. Thus, leptons which are initially above (be- 
low) the (x, z) plane drift down (up) towards this plane as 
sketched in Figure 1. A restricted class of solutions has 
the form x(t) = [x(t), 0, z(t)\. A particle "launched" from 
a position (xo, 0, 0) in the y = plane with v y = stays in 
this plane. Of course, a particle launched away from this 
plane (or one that has v y ^ 0), oscillates about the y = 
plane. 

In this work we neglect the thickness of the current layer 
in the y direction. Therefore, we use the particle orbits in 
the y = plane. We consider a set of particles entering the 
y = plane from above and below with initial locations 
given by (x ,0,0). 

Dimcnsionless variables are introduced as 

x = x/L, i = Ldot , where Ld n = \q\B /mc , (7) 

is the non-relativistic cyclotron frequency in the reference 
field B - Dimensionless electric and magnetic fields are 
naturally measured in units of _B - In particular, we let 

£=\E \/B Q . (8) 

A characteristic Lorentz factor can be defined as 
\q\B L uj L 



7o = 

We can rewrite this as 
7o w 5.9 x 10 3 



me 1 



(9) 



B n 



10- 6 G 7 V 10 13 cm 



Note that v x /c = -f (dx/dt), for example. We let be 
the number density of leptons (electrons + positrons) at 
a large distance from the X— point. An important dimen- 
sionless quantity is 



c 



47m oo 7 mc 2 



(10) 



The role of this parameter is described in §3. 

2.3. High Energy Approximation 

A particle launched at xq in the y = plane is acceler- 
ated in the z— direction by the E z field and it drifts in the 
x— direction because of the v x B force. The closer to the 
X— point a particle starts (the smaller xo is) the higher 
the energy it is accelerated to. An approximate solution 
for the particle motion in the external magnetic field B ea:t 
is possible for sufficiently small x where there is a large 
acceleration in the z— direction. 

Assuming that the particle is already highly relativistic 
we let p x (t) = m-fcsin(9) and p z (t) = mjccos(9). For 
specificity we consider the acceleration of electrons. The 
equations of motion become 

' h '- -qB y *\x) cos(8) , ^ = qB^{x) sm(9) + qE . 



We assume 2 -C 1 and discuss the conditions for this later. 
The energy equation gives mjc 2 = qE^z + const and the 
small angle approximate gives t ss z/c. 

The x— equation of motion in dimensionless variables 
becomes 
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dz 
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The relevant solution to equation (11) is 



(11) 



(12) 



where To is the usual modified Bessel function. This depen- 
dence is valid for small angles, that is, for (dx/dz) 2 <C 1, 
or [xoh(2^y z/£)/\fi£] 2 < 1. A necessary condition for 
this is that xo <C £ ■ The geometry of the orbits for highly 
relativistic motion depends only on xo and £ (in fact, on 
just £/xq) and not for example on 70. 

Figure 2 shows orbits in the external magnetic field cal- 
culated without approximation. The approximation of 
equation (12) is also shown. 

2.4. Current- Density from Single Particle Motion 

Far from the A— point, ideal magnetohydrodynamics 
applies, and the fluid particles are confined to flux surfaces 
where A z =const. Even though this is a time independent 
situation, the flux surfaces drift inward towards the neu- 
tral layer with a velocity equal to = cE x B/£> 2 . Since 
the B field is tangent to surfaces of constant A z and the 
electric field is in the z direction, the field lines drift in- 
ward towards the X— point at the same speed the particles 
do. Near the X point, where |E| > |B|, ideal magnetohy- 
drodynamics breaks down. 

We assume that after a particle inflows (from above or 
below) to the y = plane it gets launched on an orbit in 
this plane, [x(t), 0, z(t)]. The initial position of the particle 
in the y — plane is xo and some value of z. We denote 
the number flux density of particles (electrons + positrons) 
inflowing from above and below the y — plane as F{x$) 
(number per cm 2 per second). The simplest case is that 
of a uniform density of inflowing plasma where 



T = 2nooc£ 



(13) 



dt 



where is the number density (1/cm 3 ) of electron plus 
positron density at large distances and where the factor 
of two comes from the two sides of the current layer. The 
motion of a single particle starting from position x is de- 
scribed by its position [a:(£|xo), 0, z(t\x )} and its velocity 
[v x (t\x ),0,v z (t\x )]. 

The surface current-density (charge per cm per second) 
for each initial xq is then given by 

J z (x\x )dx Q = qAf(x\x Q )(v z (x\xo))dxo , (14) 

where Af(x\xo) (with units 1/cm 3 ) is the surface number 
density (of electrons and positrons) at x launched between 
xo and xo + c?xo. In equation (14), v z is the electron ve- 
locity, q is the electron charge, and the angular brackets 
indicate a time average. In the considered stationary state 
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Fig. 2. — Sample electron orbits in the y = plane for the external magnetic field By Xt = Bo x/L for the case £ = 0.5, 
7o = 10 3 , and oto = 0.1, 0.2, and 0.3. For large x, where the orbits are "looping," the drift in the +x— direction is the 
E z x By drift, and the drift in the +z direction is the gradient B drift oc B y x VB y {x, 0). The dashed line segment is the 
approximate dependence of equation (12). The positron orbits are mirror images for z — > — z. 



which is uniform in z, all averages are independent of t and 
z. 

Consider particles moving in the +x direction (the —x 
direction case follows by symmetry) , then conservation of 
particles implies that 



N(x\xq){v x (x\xo)) = const = !F{xq) 



(15) 



where v x {x\xq) is the velocity of an electron or a positron. 
Thus, 

Jz{x) = I dx J z (x\x ) , 
Jo 

= q dx T(x ) , 16 

Jo (Vx{x\x )} 

where X™ is the maximum of Xq assumed to be 0{L). The 
total current between x and x + Sx can be written as 



5x J z {x) = q dx T{xq) 
Jo 



X+Sx , (v t ) 
ax -. — r 



For a particular value of Xq a particle orbit passes be- 
tween x and x + Sx for time between t\ = t(x\xo) and 
t 2 = t(x + Sx\xq) (and possibly for time between i 3 and 
ti and ts and tg, or an odd number of intervals). This 
allows us to replace the spatial integration by a tem- 
poral integration over the particle orbit. Notice that 
(vx) — ^2$tjV x j/ "%26tj, where j = 1 or j = 1,2,3, etc. 
represents the different traversals of the interval 5x at the 
distance x, and where Stj — Sx/\v x j\. We have for example 
J dtv x = StjV X j = Sx. Therefore, we find 



f x 

Sx J z (x) = q dx F(xo) / dt v z [x 
Jo Jti 



(t),z(t)\x ] • 



(16a) 

The time integral needs to be done for all values of t for 
which the particle is between x and x + Sx. We can put 
equation (16a) into dimensionless form using equations (7) 
- (10) and J z = J z /{cBq). For the case where equation 
(13) applies this gives 



Sx J z 



= Qjdx J 



dt 



dz 
di 



Here, 



2ir ' 



(166) 



(17) 



with £ defined by equation (10), is a dimensionless measure 
of the rate of inflow of plasma to the neutral layer. 

2.5. Self-Magnetic Field 



For the considered thin current layer with surface 
current-density J{x), we can express the self- magnetic 
field of this current as 



B s y elf (x) = - c vJdx 



, MX') 



(18) 



where V is the principle value. In dimensionless variables 
this becomes 



Bf f {x) =2vJ dr 



, Jz{x') 



X — X' 



The total magnetic field is 

B y ot (x) = B c y xt (x)+B y el f(x) 



(19) 



The dimensionless form of this equation is the same but 
with hats over the magnetic fields. 

In the following we simplify the notation by dropping 
the hats and dropping the y and z subscripts on B y and 

Jz- 

3. METHODS AND RESULTS 

The self-consistent magnetic field of the current layer is 
calculated iteratively. For a relatively weak self field - suf- 
ficiently small ( compared with unity - a direct iteration 
scheme has been found to converge. This scheme may be 
described as 



Bl 



where 



Jn+l 



Tjtot 



ra+2 



(20) 



orbits 



HT 



orbits 



B'f = x + Blf , B s n elf = 2 vf dx' . (21a 

J xx 



Here, "HT" stands for the Hilbert transform, and "orbits" 
indicates the evaluation of equation (16b) using an x— grid 
of 100—200 points and the calculation of 100—200 electron 
orbits. A stretched x— grid is used with higher resolution 
at small x. The distribution of inflowing plasma !F{xq) 
is assumed uniform in x as given by equation (13). All 
cases considered here have 70 S> 1 so that the initial, non- 
relativistic or modestly relativistic conditions (7 < 2) of 
particles at their starting points x do not affect the re- 
ported results. In the iteration sequence we have found it 
advantageous to make analytical fits to the self fields B^ 
which are used in the subsequent orbit calculations. 

For larger self-fields we use the iteration (20) with equa- 
tion (21a) replaced by 



B l ° l = x + aB s n elf + (1 - a)B, 



self 



(216) 



where < a < 1 is a numerical factor, and B^ is still 
given by equation (21a). This iteration scheme is known to 
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allow the calculation of self-consistent particle ring equi- 
libria with large self-fields (see for example Larrabee et al. 
1979, 1982). 

The key dimensionless parameters which determine the 
self-consistent current layers are 



C , and 70 . 



(22) 



For the considered limit 70 » 1 and fixed £, we find that 
the current density and self-field are almost independent 
of 70. However, the distribution function of accelerated 
leptons does depend on 70 and £ as discussed below. 

3.1. Weak Self-Field 

Figure 3 shows the current density profiles J(x) for a 
weak field case C = 0.05 where the convergence of the 
direct iteration, equations (20) and (21a), is quite rapid. 
The main peak of the current density is due to the direct 
acceleration of electrons and positrons in the region where 
|E|/|B| > 1. The "shoulder" of the current density profile 
is due to the gradient B drift in the z— direction. Figure 4 
shows the final profiles of the total and self-magnetic fields. 
The self-field acts to reduce the slope of the total field B tot 
near the origin. For vanishing self-field the slope is unity. 
For the case of Figure 4, the slope is dB tot /dx\o « 0.680. 

The energy distribution or spectrum of the leptons ac- 
celerated in the current layer dn/dj is obtained from the 
calculation of the particle orbits. The particles inflowing 
to the current layer are uniformly distributed in xo and ini- 
tially they are assumed to have 7 < 2. The Lorentz factor 
of a particle as it exits the current layer (say, x > x max ) 
is therefore a function of x , "f m — 7 m (x ). As part of 
our orbit calculations we determine 7 m (xo). We find that 
7m (xo) is a monotonically decreasing function of xq. With 
dn/dxQ = const the number of particles launched between 
Xq and xq + dx , we have 



dn dn / dxo 
dj \d"f/dx \ 



(23) 



for 7 > 2, where we have dropped the m— subscript on 7. 

Figure 5 shows lepton distribution for the same case 
as Figures 3 and 4. As mentioned earlier the magnetic 
field of the current layer is almost independent of 70. For 
these values of C, and £, we have run a range of values of 
70 = 10 2 — 10 4 , and find that the distributions are roughly 
fitted by 

dn K ( 7 \ . , 

where T = k £ 2 7o and k w 12. With K w l/ln(0.56r ) 
the distribution is normalized to unity. The distribution 
(24) is not expected to be the same as that observed 
because of the influence of the distribution of 70 values 
and/or the affect of radiative losses. These affects on the 
spectrum are discussed in §3.4. 

3.2. Strong Self-Field 

For larger self-field strengths £, we use the iteration in- 
dicated by equations (20) and (21b) usually with a = 0.5. 
For ( = 0.1 we have not succeeded in getting the iterations 
to converge. We find that the iterations lead to negative 
values of B tot (x) for x > which would cause trapping of 



particles in the x— direction. Such a configuration is incon- 
sistent with our assumption that the accelerated particles 
are expelled in the x— direction. The existence of a max- 
imum of (, which measures the rate of plasma inflow to 
the neutral layer, is consistent with the finding of a max- 
imum plasma inflow rate by Burkhart et al. (1991). For 
( = 0.075, equations (20) and (21b) converge after five it- 
erations. We find that the maximum of the current-density 
J(Q) and the half-width _at half-maximum of the current- 
density Ax increases as ( increases. For this case we find 
dB tot jdx 0.412. The energy spectrum of the accelerated 
leptons is similar to equation (24). 



3.3. Scalings 
In more detail, we find the scaling relations 
J(0) w 1.53 C S/B' , 



and 



Ax : 

dB tot /dx\ x 



IA8£/VW 



(25a) 



(256) 



where B' = dB tot / dx\ x -o. The first relation can be de- 
rived from the analytic orbits of §2.3. For the external 
magnetic field B tot = x, we have from §2.4, 



lim J(x) 

x— >0 



limC 



dx 
dx/dz 

; dz' 



lim ( 



dxo 



x /i(2a 



/o(2 



1.53C£ 



(26) 



where the relation between x, xo, and z is given by equa- 
tion (12). In the magnetic field B tot — B'x, equation (26) 
is modified by the replacement £ — > £/B'. Equation (25b) 
is empirical. 

We can use equations (25a) and (25b) to obtain an ap- 
proximate constitutive equation for the current layer. Note 
that dB tot /dx\ = l+dB sel f /dx\ , and that dB sel * /dx\ = 
—2 / dx[J(0) — J(x)]/x 2 . Notice in turn that this integral 
is proportional to j(0)/Ax. From this we obtain the re- 
lation 

C w 0.205 [(B') 1/2 - (B') 3/2 1 • (27) 

This dependence is shown in Figure 6 along with the 
calculated equilibria. The left-hand part of the curve 
is dashed because we have not found the correspond- 
ing self-consistent equilibria. The maximum of ( is thus 
Cmax ~ 0.08. This maximum represents a maximum inflow 
rate to the reconnection region analogous to that found by 
Burkhart et al. (1991). 

3.4. Apparent Spectra 

The observed lepton distribution f(-y) deduced from the 
synchrotron spectrum of a radio source will in general be 
different from the distribution of accelerated particles at a 
given reconnection site (dn/d'y). For a typical radio source 
/ cx 7~ 2 - 5 . One effect is due to the distribution of r values 
owing mainly to the distribution of L for different recon- 
nection sites. Note that k = 2ir/L is the wavenumber of 
the spatial power spectrum of the B field. With the dis- 
tribution VF(ro), the average lepton distribution function 
7(7) for many reconnection sites is 



/(7) = J dT W(T ) 



dn(j\r ) 
07 



(28) 
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Fig. 3.— Iteration sequence of the dimensionless current density for a case of weak self-field where ( = 0.05. For this 
case £ — 0.5 and 70 = 10 3 . 



Fig. 4. — Dimensionless self-consistent magnetic field profiles at the fourth iteration for same case as Figure 3, C — 0.05, 
£ = 0.5, and 7o = 10 3 . The dashed B sel f C curve is the analytic fit given by 5.4a;(l + 1.9a; 2 )/(l + 0.52a; 2 + 2.1a; 4 ). Near the 
origin the slope of the total field, dB tot /dx\o ~ 0.680, is reduced from unity owing to the self-field of the current layer. 



Fig. 5. — Spectrum of lcptons accelerated in the current layer for £ = 0.05, £ = 0.05, and 70 = 10 3 and 3 x 10 3 . 
For this spectrum we have taken x max = 2.5. The dotted curve is an approximate fit for the 70 = 10 3 case given by 
dn/d'y = (K/-f) cxp(— K7/70) with k = 1/3. With K w 1/ m(0.56l7o/re), the lepton distribution function is normalized to 
unity. 



For W(T ) oc Tg" 9 , equation (28) gives / oc j- 1 ^ for 
q > 1, where the logarithmic dependence of the normal- 
ization factor K of dn/d'y is neglected. 

The observed /( 7 ) may also be steeper than dn/d'y 
owing to synchrotron and/or inverse Compton radiation 
(see, e.g., Pacholczyk 1970). In a steady state the radia- 
tive loss term in the Boltzmann equation for the leptons, 
^[7/ (7)]/^7 balances the source term dn/d'y on the right 
hand side, where 7 oc — 7 2 describes the radiative losses. 
Thus 

, . . const f°° , . dn , . 

m d7 V (29) 

This integral is easily done for the case where dn/d'y is 
given by equation (24). One finds that the f(j) is some- 
what steeper than 7 ~ 2 for 7 < T and rapidly decreasing 
for 7 > r . 

4. CONCLUSIONS 

We have calculated the self-consistent current-density 
distribution in the region of collisionless reconnection at an 
X— type magnetic neutral point in a relativistic electron- 
positron plasma. The energy distribution function of the 
accelerated leptons has also been calculated. 



We find that there are three important dimensionless 
parameters: (1) The characteristic Lorentz factor 70 = 
eBoL/m e c 2 3> 1, where Bq is a reference magnetic field 
strength and L is the length scale of the magnetic field; 
(2) The ratio of the static electric field E$ to the refer- 
ence magnetic field, £ = \E n \/B ; and (3) the strength 
parameter ( = (£ /2ir)4:Trn 00 mc 2 / B 2 , which is a measure 
of the rate of inflow of the plasma towards the reconnec- 
tion region. We find that the current-density distribution 
depends on £ and £ but only weakly on 70. The distri- 
bution function of the accelerated leptons is found to be 
approximately dn/d'y oc -y^ 1 for 7 <; Tq = k£ 2 jq where 
k w 12. The simulations of Zcnitani and Hoshino (2001) 
give a distribution oc 7 _1 for 7 < 20. The apparent dis- 
tribution function may be steeper due to the distribution 
of To values and/or the radiative losses. Self-consistent 
equilibria are found only for plasma inflow rates £ to the 
X — point less than a critical value which depends on £ . 
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Fig. 6. — Constitutive relation for relativistic electron-positron reconnection from equation (27). The solid circles are 
calculated equilibria for £ = 0.5 while the open circles are for £ = 0.2. No equilibria have been found along the dashed 
part of the curve. 
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